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Abstract 

The transition from the quantum Hall state to the insulator is considered for non-interacting electrons in a two- 
dimensional disordered lattice model with perpendicular magnetic field. Using correlated random disorder poten- 
tials the floating up of the critical states can be observed in a similar way as in the continuum model. Thus, the 
peculiar behaviour of the lattice models reported previously originates in the special choice of uncorrelated random 
disorder potentials. 
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1. Introduction 

Instead of floating up in energy, as has been pro- 
posed [pi 0] and shown to exist in continuum 
models [§], a peculiar annihilation of the cur- 
rent carrying states takes place in lattice models 
when the quantum Hall to insulator transition is 
approached. According to recent investigations 
© 1' H' El' mcreasm g tne disorder causes states 
associated with a negative Chern number to start 
moving outwards from the band centre and to 
destroy states with corresponding positive Chern 
numbers at the band edges. This behaviour is a 
special feature of lattice models which implicitly 
contain a periodic potential, but is not possible in 
the continuum model without a periodic potential 
usually considered to explain the levitation of the 
current carrying states. 



The transition from the quantum Hall phase to 
the insulator is strongly affected by the way the 
critical states disappear. While in the continuum 
model the crossover takes place when the last cur- 
rent carrying state floats across the Fermi level, 
a direct transition to the insulator is possible in 
the lattice model also from higher Hall plateaus. 
The latter possibility is not contained in the global 
phase diagram for the quantum Hall effect [ |l0| ] 
since it was based on the levitation scenario [jj],gp. 

Up to now the results of experimental investi- 
gations are not conclusive [fljl [12, 13, 14]. Both, 
direct transitions from higher Hall plateaus to the 
insulator have been reported as well as experi- 
ments that are in accord with the global phase 
diagram. Recently, it has been pointed out that 
the low field quantum Hall to insulator transition 
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might be unaccessible for the present time [15], 
because the available temperature range and sam- 
ple sizes refuse to enter the scaling regime. How- 
ever, a QH to insulator transition as a function of 
disorder strength at higher magnetic fields is more 
likely within experimental reach. 
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In the present paper we consider the possible 
origin of both the diverging experimental and the- 
oretical results reported so far. In particular, the 
fate of the current carrying states is studied when 
the disorder is increased. We also investigate, 
how the QH to insulator transition depends on the 
correlation length of the disorder potentials. 

2. Model with correlated disorder potentials 

A two-dimensional tight binding model describ- 
ing non-interacting electrons moving in a perpen- 
dicular magnetic field and random disorder poten- 
tials associated with the lattice sites k is given by 
the Hamiltonian 

H = Y,w k \k)(k\+ Vki\k)(l\. (1) 

k <k^l> 

The transfer terms Vki connect only neighbour- 
ing sites and contain the magnetic field B = 
ah/ea~ 2 which is taken to be commensurate 
with the lattice. It is expressed as the number 
a = p/q of flux quanta h/e per plaquette of 
area a 2 and enters via the complex phases, Vki ~ 
exp(i27rak x /a). We take p/q = 1/8 so that the 
tight binding band splits into 8 sub-bands. The 
correlated disorder potentials Wk are generated it- 
eratively site by site using a generalisation of the 
Gibbs representation of Markov random chains 
[0]. 

The Gibbs distribution for the configuration of 
random variables fi is 



P(ji) = Z- 1 exp(-CG( M )), 



(2) 



where Q(p) = — J2<ki> a k(^) a i(l J ') is the Gibbs 
potential and Z the partition function. Here, 
crfc(/i) = /ifc G [—1/2,1/2] is the value of the 
random variable on site k. For a given configu- 
ration /i° the conditional probability that the new 
variable on the new site n takes a value between 
-1/2 and A < 1/2 is 



P( M n < A I ft ) 



-l/2C£ n 



2sinh(l/2CI7 n ) ' 



(3) 




Fig. 1. Gray-scale plot for a 128x128 array of correlated 
disorder potentials with correlation parameter C = 1.3. 



where S n = J2ieN a i> anc ^ N n is the set of 
neighbours of n already determined. The corre- 
lation parameter C tunes the strength of the po- 
tential correlations. C = gives the uncorrelated 
case. 

The correlated random numbers are then multi- 
plied by the disorder strength W which results in 
the set of disorder potentials wu used in Eqn. ([]]). 
An example of correlated potentials generated by 
this method with C = 1.3 is shown in Fig. [j] 
where a gray-scale plot of a 128 x 128 array is dis- 
played. White regions correspond to wt ~ +W 
and black areas to wu ~ — W with the gray spots 
for values in between. 

To further characterise the disorder poten- 
tials, we have calculated the correlation function 
K{p) = (wkWi) p=\r k -n\ averaged over all pairs 
of sites which are a given distance |r& — 77 1 apart, 



K(p) 



P) 



Ef^) 1/2 fef^ 1/2 ' 



(4) 



For C = 0.5,0.7,1.0,1.3,1.7, and 2.0, we find 
an exponential form, K(p) ~ exp(— — ri\/r]). 
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Fig. 2. The exponential decay of the correlations in the dis- 
order potentials for C = 0.5 (+), 0.7 (x), 1.0 (*), 1.3 (□), 
1.7(B), and 2.0 (O). 



The spatial decay of the correlations is governed 
by the correlation length 77 which in the range 
0.7 < C < 2.5 can approximately be fitted by 
77(C) Ki exp(2C). The results for the correla- 
tion function K(p) and the correlation length 77 
describing the decay of the potential correlations 
are shown in Fig. ||. 

3. Hall conductivity and localisation length 

The zero temperature Hall conductivity 



&xy(E, W) 
L 



e 2 

lim lim TW (5) 

e->0+ L^oo h U 



J>(G± - G- l )x t y-2Y,e 2 G±yGj l xA 

i i,j 



and the localisation length 



lim _lnTr|Gi, L |^ (6) 

L— >oo JjLi 



are calculated numerically for large disordered 
two-dimensional systems of width M and length 
L using a recursive Green function method de- 
veloped previously [17, [l^]. Depending on the 
correlation length 77 the width is varied between 
M/a = 32 and M/a = 160. The length of the 
system necessary for convergence is typically be- 
tween 10 5 and 10 7 in units of the lattice constant 
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Fig. 3. Normalised localisation length Xm(E, W)/M as a 
function of energy E for a system width M/a = 48. The 
correlation parameter is C = 0.2. The disorder strength is 
W/V = 0.5 (O), 1.5 (+), 2.0 (□), 2.2 (x), and 2.5 (A). 



a. In calculating the Hall conductivity one has 
strictly to observe the correct order of limits. 



4. Results and discussion 

The Hall conductivity was calculated near filling 
factor v = 2 as a function of increasing disor- 
der W. For correlation parameters in the range 
0.0 < C < 1.0 a continuous decay from a xy = 
2e 2 /h —>■ a xy = is observed. Although a direct 
transition is seen in the Hall conductivity for these 
correlation parameters, a closer look at the cor- 
responding localisation length indicates a totally 
different behaviour of the current carrying states. 
For C = 0.2, which corresponds to a correlation 
length smaller than the magnetic length, no float- 
ing up of the critical state of the lowest Landau 
band can be observed. In fact, as in the uncorre- 
cted case it slightly moves downward in energy 
with increasing disorder strength until it gets an- 
nihilated by the corresponding state with negative 
Chern number moving faster downward from the 
band centre. This is shown for a system of width 
M/a = 48 in Fig. [| where the normalised local- 
isation length Xm/M is plotted versus energy in 
the filling factor range < v < 1.5 and several 
disorder strength W. The transition to the insula- 
tor takes place at W c /V 2.5 where the Chern 
and anti-Chern state touch each other. 
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Fig.4. Normalised localisation length \m(E, W) /M versus 
strength of the disorder potentials W. C = 1.0 and M = 
48 a. The energy parameter is E/V = —1.2 (O), —1.3 
(+), -1.4 (□), -1.5 (x), -1.6 (A), -1.7 (*), and -1.8 (O 
dashed line). 



A completely opposite behaviour is obtained 
for a larger correlation parameter, C = 1.0. This 
corresponds to a correlation decay length of 1.9 a 
which is almost twice as large as the magnetic 
length l B /a = (27ra) -1 . While for C = 0.2 
the last current carrying state disappeared at a dis- 
order strength of about W/V 2.5 two critical 
states are still observed at W/V k 3.5 and C = 
1.0. This is shown in Fig. ^ where the disorder 
dependence of the localization length for energies 
corresponding to filling factors 1.5 < v < 2.5 is 
displayed. Fig. || clearly reveals that there still 
exist two distinct critical points which move to 
higher disorder values when the filling factor is 
increased. In other words, the current carrying 
states of the two lowest Landau bands float up in 
energy and cross the Fermi energy one by one. 
Within this energy (filling factor) range the shift 
in disorder is weaker for the critical state of the 
lowest Landau band while it is stronger for the 
next critical state so that the two states get closer 
with increasing filling factor. 

Since the two critical points are close in dis- 
order the calculated Hall conductivity even for a 
system of width M/a = 160 does not reflect the 
situation of two seperate transitions. The reason 
for this behaviour is that the divergence of the lo- 
calization length of the lowest and that of the next 



current carrying state overlap because for a given 
Fermi energy the two critical disorders are close 
(see abscissa in Fig |j) so that the Hall steps can- 
not be resolved. This overlap will eventually dis- 
appear in the thermodynamic limit and similarly 
will the Hall conductivity exhibit the missing Hall 
plateau at a xy = e 2 /h. 

In conclusion, a lattice model with sufficient 
correlation in the disorder potentials shows sim- 
ilar behaviour as the continuum model. Hence, a 
lattice model can perfectly be in agreement with 
the proposed levitation scenario ||] and the pre- 
dictions of the global phase diagram [|lG|]. 

References 



[1] 

[2] 
[3] 



R. B. Laughlin, Phys. Rev. Lett 52 (1984) 2304. 
D. E. Khmelnitskii, Phys. Lett. 106A (1984) 182. 
T. Ando, Journal of the Physical Society of Japan 53 
(1984) 3126. 

[4] V. Kagalovsky, B. Horovitz, and Y. Avishai, Phys. 

Rev. B 52(1995)R17044. 
[5] D. Liu, X. Xie, and Q. Niu, Phys. Rev. Lett. 76 (1996) 

975. 

[6] D. N. Sheng and Z. Y. Weng, Phys. Rev. Lett. 80 
(1998) 580. 

[7] H. Potempa, A. Baker, and L. Schweitzer, Physica B 

256-258 (1998) 591. 
[8] K. Yang and R. N. Bhatt, Phys. Rev. B 59 (1999) 

8144. 

Y. Hatsugai, K. Ishibashi, and Y. Morita, Phys. Rev. 
Lett. 83 (1999) 2246. 

S. Kivelson, D.-H. Lee, and S.-C. Zhang, Phys. Rev. 
B 46 (1992) 2223. 

A. A. Shashkin, G. V. Kravchenko, and V. T. Dolgo- 
polov, JETP Lett. 58 (1993) 220. 
[12] I. Glozman, C. E. Johnson, and H. W. Jiang, Phys. 
Rev. Lett. 74 (1995) 594. 

S. V. Kravchenko, W. Mason, J. E. Furneaux, and 
V. M. Pudalov, Phys. Rev. Lett. 75 (1995) 910. 



[9] 
[10] 
[11] 



[13] 



[14] M. Hilke era/. (1999) cond-mat/9906212 



[15] B. Huckestein, Phys. Rev. Lett. 84 (2000) 3141. 
[16] R. Johnston and L. Schweitzer, in High Magnetic 
Fields in Semiconductor Physics: Proc. of the Int. 
Conference, Wiirzburg 1986, ed. by G. Landwehr 
(Springer Series in Solid-State Sciences 71, Springer, 
Berlin, 1987), pp. 71-74. 
[17] A. MacKinnon, Z. Phys. B 59 (1985) 385. 

L. Schweitzer, B. Kramer, and A. MacKinnon, Z. 
Phys. B 59 (1985) 379. 



[18] 



4 



